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Abstract 

A novel strategy to handle divergences typical of perturbative calculations is implemented for the 
Nambu- Jona-Lasinio model and its phenomenological consequences investigated. The central idea 
of the method is to avoid the critical step involved in the regularization process, namely the explicit 
evaluation of divergent integrals. This goal is achieved by assuming a regularization distribution 
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in an implicit way and making use, in intermediary steps, only of very general properties of such 
regularization. The finite parts are separated of the divergent ones and integrated free from effects 
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ly-j , of the regularization. The divergent parts are organized in terms of standard objects which are 



independent of the (arbitrary) momenta running in internal lines of loop graphs. Through the 
analysis of symmetry relations, a set of properties for the divergent objects are identified, which we 
denominate consistency relations, reducing the number of divergent objects to only a few ones. The 
calculational strategy eliminates unphysical dependencies of the arbitrary choices for the routing 
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of internal momenta, leading to ambiguity- free, and symmetry-preserving physical amplitudes. We 
show that the imposition of scale properties for the basic divergent objects leads to a critical 
condition for the constituent quark mass such that the remaining arbitrariness is removed. The 
model become predictive in the sense that its phenomenological consequences do not depend on 
possible choices made in intermediary steps. Numerical results are obtained for physical quantities 
at the one-loop level for the pion and sigma masses and pion-quark and sigma-quark coupling 
constants. 

PACS numbers: 12.39.-x,12.38.Bx,11.30.Rd 
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I. INTRODUCTION 



It is largely accepted that many of the essential features of chiral symmetry in quantum 
chromodynamics (QCD) are captured by the Nambu-Jona-Lasinio (NJL) model [lj], a simple 
relativistic quantum field theory with (nonrenormalizable) four-fermion interactions. In the 
limit of exact chiral symmetry the fermions are massless and the interaction Lagrangian 
density of the model, in its simplest version, contains the chirally symmetric sum of scalar 
and pseudo-scalar four-fermion interactions. Since the first works, using the model with 
quark degrees of freedom, in the earlier days of QCD 

aa 

, the model has been extensively 
used to describe low energy hadronic observables, like hadronic masses, correlation and 
structure functions in vacuum and at finite densities and temperatures - for a complete list 



of references see the reviews in Refs. 



Qflflflfly, 



id, 



One of the reasons for the widespread use of the model is that it realizes the dynami- 
cal breaking of chiral symmetry already at the one-loop (mean field) approximation. The 
predictions of the model, however, are intimately compromising with the specific strategy 
adopted to handle the ultraviolet divergences given the nonrenormalizable nature of the 
model. As a consequence, the specification of a procedure for handling divergent amplitudes 
is a necessary and essential first step to be performed before extracting physical predictions 
as must be made in any relativistic quantum field theory. In a renormalizable theory this 
is done by specifying a regularization procedure by which the divergences are isolated and 
eliminated through a convenient reparametrization of the theory, removing, in this way, 
any parameters introduced in the regularization process. Therefore, although essential and 
necessary, the regularization process plays a secondary role and seems to be a disposable 
intermediate step, in the sense that it is not meant to modify the physical content of the 
theory. However, there is a distinctive and nontrivial aspect with the regularization of the 
ultraviolet divergences in the NJL model in view of its nonrenormalizability. Unlike with 
renormalizable models, as an increasing number of loops is considered the reparametrization 
of the model can be made only at the cost of adding an increasing number of terms with 
extra coupling constants to the original Lagrangian in order to render physical amplitudes 
independent of the regularization procedure. In principle, there is no problem with a theory 
having an infinite number of coupling constants when using it as an effective field theory, 
as explained by Weinberg [12|. However, practitioners of the NJL model have followed 
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the attitude of using it as a regularization-dependent model, considering the regularization 
procedure part of the definition of the model. Within such an attitude a large body of 
interesting and valuable work has been, and continues to be done using the model. 

The regularization of divergent amplitudes is a delicate process due to the arbitrariness 
in the manipulation of improper integrals which can be converted into ambiguities when 
results become dependent on the choices involved. There are ambiguities associated with 
the arbitrary routing of the momenta in internal lines of divergent loop amplitudes, which 
invariably lead to the violation of space-time homogeneity. There are also ambiguities asso- 
ciated with the choice of the common scale for the divergent and finite parts of amplitudes 
that may lead to the breaking of scale invariance. In general, different sorts of ambiguities 
have the potential of leading to violations of symmetry relations of global and local gauge 
symmetries. The most commonly used regularization procedures for the NJL model such as 
three- and four-momentum cutoff, Pauli-Villars and proper-time lead to one or more of such 
symmetry violations. Dimensional regularization (DR), although not much used within the 
NJL model, in general, leads to amplitudes free from ambiguities and symmetry preserv- 
ing. However, it has problems at high densities and temperatures, when chiral symmetry 
is restored. This is due to the fact that the quadratic divergence which appears in almost 
all one- loop amplitudes must be assumed as zero in the zero-mass limit. Practice with the 
NJL model has shown that depending on the problem studied, one regularization scheme 
seems to be more appropriate than another because of the problems just mentioned. For 
example, when working with correlation functions, in general dispersion relations are not 
automatically fulfilled in cutoff and proper-time regularizations, in contrast to Pauli-Villars 
regularization. On the other hand, while causality is preserved with Pauli-Villars regulariza- 
tion, unitarity is violated at high enough energies, although it is preserved with proper-time 
regularization. These sorts of problems are well known and arguments have been put for- 
ward and tricks invented to deal with such problems - for a discussion on these issues, see 



for example Refs. 13j, 



151 ] . Obviously, this situation is unsatisfactory since one would like 



that the regularization scheme play a secondary role in the process of making predictions 
with the model. 



The difficulties pointed out above lead Willey [1G] and Gherghetta [17| to conclude that 
there is no way to make consistent physical predictions with the NJL model using tradi- 
tional regularization techniques. However, this question in the context of the gauged NJL 
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model was considered in a later work by Battistel and Nemes 18| using a novel strategy 



to handle divergent amplitudes 



19j. The referred investigation revealed that NJL ampli- 



tudes ambiguity-free and symmetry-preserving can be obtained, and as such making the 
NJL model predictive. The central idea of the method is to avoid the critical step involved 
in the regularization process, namely the explicit evaluation of divergent integrals. This goal 
is achieved by assuming a regularization distribution in an implicit way and making use, in 
intermediary steps, only of very general properties of such regularization. The finite parts 
are separated of the divergent ones and integrated free from effects of the regularization in a 
completely similar way as made in the treatment of renormalizable theories. The divergent 
parts are organized in terms of standard objects which are independent of the (arbitrary) 
momenta running in internal lines of loop graphs. Through the analysis of symmetry re- 
lations, a set of properties for the divergent objects are identified, which we denominate 
consistency relations (CR's), reducing the number of divergent objects to only a few ones. 
The remaining objects never really need to be evaluated. In renormalizable theories they 
are eliminated by the counter-terms. In a non-renormalizable model, such as in the NJL, 
the basic divergences are fixed by fitting observables, as we will see along this contribution. 

Having in mind this perspective, in the present work we extend the original discussion of 
Ref. [l8[ by presenting a complete and unified discussion on the symmetry relations involving 
Green's functions, including tensor operators. Some of the relations can be derived using 
the methods of current algebra, in particular when using the conservation of the fermionic 
vector current and the proportionality of the divergence of the axial vector current to the 
pseudoscalar current. Here we also discuss relations of Green's functions of tensor operators. 
Such relations cannot be obtained through current algebra methods, because the divergence 
of the fermionic tensor current cannot be written in terms of other fermionic currents. This 
is the case, for example, of the relations involving tensor-tensor two-point amplitudes. In a 
first step of our investigation we will show that it is possible to obtain physical amplitudes 
preserving the symmetries and automatically free from ambiguities associated with the ar- 
bitrariness in the routing of momenta in the internal lines of loops. In a second step we 
will show that the undefined quantities associated with the divergent objects can be fixed 
phenomenologically, without to the recourse of calculating any divergent integral, leading in 
this way to a regularization independent parametrization of the model. The model, within 
this formulation, becomes predictive in the sense that all arbitrariness are removed from the 
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physical amplitudes. The model works like a "renormalizable" theory at the one-loop level. 

The results obtained in the present paper are new and extend the applicability of the 
NJL model in way that it becomes independent of a particular regularization scheme, since 
no explicit regularization is actually used. It is new because all symmetry constraints on 
general Green's functions, including tensor ones, are preserved and as such no problems 
with causality and unitarity can arise. The present method has also been applied in the 
context of CPT breaking in models with Chern-Simons interactions [20], to the neutral 
electromagnetic pion-decay where the AVV triangle anomaly phenomenon is discussed, tri- 
angle anomalies 22j and tensor densities 23]. In particular, it was shown that the adopted 
strategy furnish the expected anomalous term and the ambiguities again play no relevant 
role 21(. One important aspect of the method presented here is that, since no explicit 
regulator is used, a result obtained within a given traditional regularization method can be 
immediately reobtained by explicit evaluation of the implicitly regulated expressions. With 
this, the use of the consistency relations allow to identify the reasons why commonly used 
regularization schemes lead to symmetry violations. Invariably, the reasons are that not all 
the consistency relations are satisfied simultaneously within the traditional regularization 
methods. 

The plan of this paper is the following. In Section [Til we present the Lagrangian density 
of the NJL model used in this paper, and discuss general "consistency constraints" that the 
one-loop Green's functions must satisfy in order not to violate symmetry relations. Next, 
in Section II II I we discuss a calculational scheme that preserves the general relationships 
among the Greens functions obtained in Section [TTJ The calculational scheme isolates the 
purely divergent terms, which will disappear because of symmetry consideration or will 
be fitted to observables, while the finite parts are integrated without any regularization. In 
section HV] we discuss the ambiguities associated with the choices for the momentum routing 
in the internal lines of loops and show that the methods used to isolate the divergent parts 
respects all the general relations among Green's functions. This aspect of the regularization 
is central to the paper and is highly nontrivial since, because of the ultraviolet divergences, 
these relations can very easily be violated when not being careful with the explicit evaluation 
of the integrals within a particular regularization scheme. In Section [V] we identify the 
general properties that the divergent ambiguous quantities must satisfy in order to guarantee 
the preservation of fundamental symmetries. For that, Ward identities and other general 
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constraints imposed by Furry's theorem will be invoked. The phenomenology and numerical 
results are presented in Section [VTl Here only the traditional observables, like pion and sigma 
masses and coupling constants are calculated to show that the traditional phenomenology 
is obtained in a straightforward way Our Conclusions and Perspectives for future work are 
presented in Section IVHl 



II. MODEL LAGRANGIAN AND CONSISTENCY CONSTRAINTS 

In this paper we restrict the discussion to the simplest SU(2) version of the NJL model 
that incorporates the light-quark u and d flavors only The SU(3) case will be considered 
elsewhere [24J. The Lagrangian density is given by 

£ = -m )ip + G S (tptp) 2 + (ipri^ifj) 2 , (1) 

where ip is the two-flavor, three-color quark field operator and m is the diagonal current 
quark mass matrix. To simplify the discussion, we take equal u and d quark masses. The 
nonperturbative quark propagator S{p) is given in terms of the self-energy S(p) as 

s-\ V )=^-nv)- (2) 

In the mean field approximation, the self-energy is momentum independent £Qo) = M, with 
M satisfying a gap equation [l| 



M = m -2G s N f (W), (3) 

where Nf = 2 is the number of flavors, and (ipip) is the one-flavor, Lorentz scalar one-point 
function (the quark condensate) given by 

-Tr[S(k)} = —AN c t / — --- — , (4) 



(2tt) 4 l v /j c J (27r) 4 P-M 2: 

where N c = 3 is the number of colors. 

In general, phenomenological predictions for meson masses and correlation functions re- 
quire the evaluation of purely fermionic n-point Green's functions. In spite the fact the 
model is nonrenormalizable, the Green's functions obey well-defined relations among them. 
Such relations are the manifestation of the symmetries of the underlying Lagrangian defining 
the model. Therefore, in any attempt of describing a specific phenomenology it is crucial 



that the evaluation of physical amplitudes preserve such symmetry relations. If it turns 
that such symmetry relations are not preserved by the calculation, the predictions cannot 
be characterized as consequences of the underlying symmetries supposed relevant for the 
specific phenomenology and which were the main motivation for using a schematic model 
like the NJL model. 

A generic one-loop n-point Green's function can be defined as 



T rir2 "- r " (k u k 2 , ■ ■ • , k n ) = / -— -Ti [T a S (k + h) T 2 S (k + k 2 ) ■ ■ ■ T n S (k + k n )\ , (5) 



where the trace is over Dirac indices only, and the I\ represent one or more of the matrices 
(1, 7 M , 75, 7^75, a^ u ) to which we attribute the labels (S, V, P, A, T), respectively. The fermion 
propagators are given by Eq. (j2j) and the k n are arbitrary routing momenta in the internal 
lines, and are related to the external momenta. A physical amplitude can depend only on 
differences of k n , any dependence on sums of k n is unphysical. We note that the highest 
superficial degree of divergence is cubic, and occurs for the amplitude with n — 1. For n > 4, 
the amplitudes are finite. In particular, in the NJL model the scalar one-point function is 
relevant for the gap equation, the two-point functions appear in the bound-state equations 
for mesons, the amplitudes with n = 3 describe meson-decays and the four-point functions 
are relevant for meson-meson scattering. 

At the one loop level, there appear only two divergences, a quadratic and the logarithmic. 
The cubic divergence is absent in the one-point functions either because of the trace or 
because the integral is identically zero - for the same reason there is no linear divergence 
in two- and three-point functions. The standard procedure to deal with the divergences is 
to cutoff the momentum integrals at some momentum A. With this, the model has two 
unknowns, the coupling Gs and A. These can be fitted by using the values of the quark 
condensate {iptp), related to the quadratically divergent scalar one-point function T s , and 
the pion decay constant f n , related to the logarithmically divergent axial- vector pseudo- 
scalar two-point function T AP . Since these two types of divergences are the only ones that 
appear in all other Green's functions of the model, all divergences can be absorbed by the 
physical quantities (ipip) and f n . In a certain sense, this is a type of renormalization. 

There are two aspects we would like to note with respect to this "renormalization". First, 
one is explicitly using relations among divergent Green's functions having different numbers 
of points. Second, given the strict nonrenormalizability of the model, divergent amplitudes 
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are related to physical quantities through a regularization function or, in last instance, 
through the adjustment of regularization parameters which are interpreted as cutoffs in the 
momentum integration. It is well known that these two facts lead in general to symmetry 
violations. Moreover, the effects of the modification introduced in the regularization process 
remain present even in the finite parts of the amplitudes. In this respect, the regularization 
and renormalization procedure for nonrenormalizable and renormalizable models are treated 
in completely different manner. This discussion emphasizes the difficulties in satisfying 
kinematical constraints and symmetry relations involving one or more divergent amplitudes. 
The method we present here is based on the following strategy: (1) first, all constraints and 
symmetry relations are imposed without evaluating any divergent integral, (2) divergences 
are isolated into quantities that are independent of arbitrary routing momenta, and (3) 
amplitudes respecting all symmetry relations and free from ambiguities are obtained. All 
this is achieved without compromising with a particular regularization scheme. 

Let us start considering the general "consistency constraints" that Green's functions must 
satisfy in order not to violate symmetry relations. One very general and powerful way to 
generate relations among Green's functions is to identify identities, at the integrand level, 
resulting from the contraction of the Lorentz vector indices of a vertex operator with an 
external momentum. Although this method is entirely equivalent to the current algebra 
technique for some types of amplitudes, it can be applied also to tensor currents - for which 
the methods of current algebra are not applicable. As an example consider the identity 

(h - k 2 y i^S (k + h) lv S (k + k 2 )} = llx S (k + k 2 ) - lfl S (k + h) , (6) 

where S(k) is the mean field quark operator given by Eq. (j2j). This identity follows trivially 
from the algebra of the Dirac 7 matrices. Tracing both sides and integrating in momentum, 
a genuine relation among Green's functions of the model is obtained 

(h - k 2 fT v J (k u k 2 ) = Tp(ki) - Tj(k 2 ), (7) 

where and TTf are respectively the vector one-point function (Ti = 7^) and the two- 
point vector- vector amplitudes (T 1 = 7^,r 2 = j u ) - see Eq. fl5]). In the same way, one 
obtains 

(k 1 -k 2 rT^ v (k 1 ,k 2 )=T^k 1 )~T^(k 2 ). (8) 
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Similarly, one can contract a Lorentz axial-vector two-point density with an external mo- 
mentum as 

(h-hr^Sik + k^^Sik + h)] = 2M[^ 5 S(k + k 1 ) l5 S(k + k 2 )\ 

+ lv S{k + k 2 ) + ivi 5 S(k + fe 1 )75, (9) 

where we have used the anticommutation of 75 and the 7^ matrices. Again, taking the traces 
and integrating on both sides, we obtain 

(h - k 2 yr AA (h, k 2 ) = 2mt ap (k u k 2 ) + r^(k 2 ) - TZ{k x ), (10) 

where the indices A and P stand for axial and pseudoscalar corresponding respectively to 
T = 7^75 and T = 75, as defined in the paragraph following Eq. (jSJ). Following the procedure 
described above, we can get also the relations 

T S {k 2 ) — T s (ki), (11) 

-2MT PP {h, k 2 ) + T s (k 2 ) - T s {h), (12) 

-2MT PV {h, k 2 ) + T A {k 2 ) - T A (h), (13) 

T A (k 2 ) — T A (ki). (14) 

The expressions (fTll) - ffT4"l) are nothing more than the relations that follow from current 
algebra methods when using the conservation of the fermionic vector current and the pro- 
portionality of the divergence of the axial vector current with the pseudoscalar current. The 
expressions involving the one-point functions correspond to the current commutator terms. 
Additional relations can be identified at the trace level 

T$f (fci, k 2 ) = (h - k 2 ) x (T SV Y (h, k 2 ) , (15) 

Tf (k u k 2 ) = ~ (h - k 2 )^ [T ss (k u k 2 ) + T pp (k l7 k 2 )] , (16) 
T™ (h, k 2 ) - T AA (h, k 2 ) = g, u [T ss (h, k 2 ) + T pp (hi, k 2 )} . (17) 

One of the advantages of the method used above is that it can be used to obtain relations 
for Green's functions involving tensor operators. Such relations cannot be obtained through 
current algebra methods simply because the divergence of the fermionic tensor current cannot 
be written in terms of other fermionic currents. This is the case of relations involving tensor- 



(h-k 2 yT: s (k u k 2 ) = 
(h-k 2 yT AP (k u k 2 ) = 

(h-hyT^ik^h) = 
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tensor two-point amplitudes, like for example the following one 



(h - hfT^ (fci, k 2 ) = -g au \T; (k 2 ) - T v p (h) + {hi - k 2 ) p T bb (h, k 2 

[Tl (k 2 ) - T v a (h) + (h - k 2 ) a T ss (h, k 2 )] 
+ (h - k 2 ) a T AA (h, k 2 ) - (h - k 2 ) p T AA (k u k 2 ) 
-2MTZ(h,k 2 ). (18) 

Given the fact that three other Lorentz indexes are left uncontracted in the TT two-point 
function, it is also possible to establish constraints on successive contractions with the ex- 
ternal momenta. It is immediate to note that when these contractions involve both Lorentz 
indexes of a tensor operator, the result must vanish identically 

(*i - k 2 ) a (h - k 2 fT T J ap (h, k 2 ) = 0. (19) 

This is due to the fact that TjJ^p ~ a ap , where a a p = — i/2 [y a , 7/3], and therefore 
(ki — k 2 ) a (ki — k^ 13 <T a/ 3 = 0. This property imposes additional constraints on the consis- 
tent evaluation of TT two-point function. Even if such requirements seem to be obvious at 
this point, the divergent character of the integrals defining the amplitude makes satisfaction 
of this property far from being trivial. 

Following strictly the same procedure outlined above, relations involving the remaining 
tensor two-point functions can be established as 

(h-k 2 rT^ v a (k u k 2 ) = 0, (20) 
(h - k 2 fT T J (h, k 2 ) = -2MT^ P (h, k 2 ) , (21) 

(*i - (k u k 2 ) = 2MT T J (k ± , k 2 ) , (22) 

(h-hYT^Jkuh) = ~2MT^ A (k 1 ,k 2 ) = -ts uaXi :(k 1 -k 2 ) x (T sv Y(k 1 ,k 2 ), (23) 

(h - k 2 YT T J (k u k 2 ) = ±.{(k 1 - k 2 ) 2 T v u s (hi, k 2 ) - (h - k 2 ) v [T s (k 2 ) - T s (h)} } 

= - l ~e auK (h - k 2 ) a (T AV ) M (h, k 2 ) , (24) 

(h - k 2 YT T J a (h, k 2 ) = ±- (h - k 2 f [T™ (h, k 2 ) - T AA (h, k 2 )] 

- (ki - k 2 ) v T PA (ki, k 2 ) . (25) 
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In deriving these results, we have used the relations given in Eqs. (|T5|) - (|T7|) . 

At this point it is very important to note that the results obtained above for the con- 
traction of amplitudes with the external momentum are of general validity. Only algebraic 
manipulations have been made, no single divergent integral has been evaluated and no 
changes of variables under integration have been made. Thus, the results obtained are not 
compromised with any type of regularization. The important question one has to face when 
evaluating the integrals defining the different amplitudes, is how to give a meaning to the 
divergent integrals without violating these relations. This will be discussed in the next 
section. 

III. CALCULATIONAL SCHEME FOR HANDLING DIVERGENT INTEGRALS 

In the preceding section we have considered the Green's functions for which one needs to 
construct a calculational scheme that preserves the general relationships among them, which 
we denominated consistency constraints. As a matter of consistency, all the constraints must 
be fulfilled without specifying special choices of the loop momenta. Divergences will appear 
in the form of improper Feynman integrals, but only a small number of them need to be 
evaluated since all the amplitudes are combinations of a few Feynman integrals. 

The traditional way to handle a divergent Feynman integral is to adopt an explicit regu- 
larization. Invariably, this amounts in modifying the original integral in a way the integration 
becomes well defined. In the context of NJL models, such modifications are commonly made 
by introducing in the integrand a distribution in the loop momentum to render the integral 
convergent. In doing so, the results of the integrals become a function of the parameters 
of the regulating distribution. In perturbatively renormalizable theories, one tries to isolate 
the purely divergent terms in order to specify the adequate counterterms for the renormal- 
ization. The parts which are independent of the regulating parameters are identified as the 
finite parts and carry the physical content of the amplitude. This means that the functions 
of the physical momenta are not affected by the regularization in the limit the regularization 
is removed. The procedure we adopt to use for the NJL model follows closely this general 
strategy, the purely divergent parts of the amplitudes will be fitted to observables, while the 
finite parts are integrated without any regularization. 

In the first step one assumes an unspecified regularization distribution G\ f (k 2 , A?) de- 



ll 



pendent on one or more "cutoff" parameters Aj, such that the original divergent integral is 
replaced by a finite one as [19] 



/ Wf m - 1 W? m G * {k - A?) s I, w? m - (26) 



l A k r Ml r a^i 

(2n) 

The generic distribution G(k 2 , A 2 ), in addition of having the obvious property of turning the 
original integral convergent, must depend only on k 2 due to Lorentz invariance and must 
have the limit 

lim G Ai {k 2 ,A 2 )=l. (27) 

A?-*oo 

which allows us to connect the regularized integral with the original one. Having assumed 
the existence of such a regularization distribution, one manipulates the integrand of the 
divergent integral in a way to isolate all the divergences in momentum-independent integrals. 
This goal can be achieved by using the identity 

where the ki is (in principle) an arbitrary routing momentum of an internal line in a loop, and 
M is the fermion mass running in the loop. The value of the integer N is the smallest integer 
that makes integrals involving this last term finite when removing the regulating function. 
In view of Eq. (127]) . the corresponding integration can be performed without restrictions and 
will be free from the specific effects of an eventual regularization. No additional assumptions 
are made with respect to the remaining divergent terms. 

When evaluating the amplitudes considered in the previous section, after taking the 
appropriate traces, it is not difficult to convince ourselves that only five divergent Feynman 
integrals will appear, namely 

^ h) -J (27T)* [(* + kl f - M 2 ] ' { J) 

it . th. Tt^\ _ f ^k (1; k^; k^k u ) 

C 2 ' 2 ' 2 ] ~ J (2^) 4 [(* + h) 2 - M 2 ] [(k + k 2 f - M 2 ] • m 
The divergent parts of these integrals can be rewritten in terms of five divergent quantities 
that we denote by D^^, A^ u , V M ^, h og and I qua d- Specifically, the integrals I\ and that 
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appear in the one-point amplitudes can be written as 



h (h) = [l q ua d (M 2 )] + k a x k{ [A af3 ] , 

h» (h) = -hfi [l qU a d (M 2 )} - fcf [V^] - ^kfk" [D a(3fMU ] 

The integrals I 2 , l 2ll and I 2fMU related to two-point amplitudes can be written as 

[/^ (M 2 )] - *(4vr)- 2 [Z {(h - k 2 f ; M 2 )] 
1 

2 1 



h (h, k 2 ) 
hp (h,k 2 ) 

huv (h,k 2 ) 



-hh + hnA^-hh + k^ih), 



2 

\2 



^ [V„„] - ^ (h - hf [A^j 

+ g ^2 ^2 + ^2 + ^1 ^l) f— U/3/J!/] 

+^ + hukg + kiuk() [A^] 

+ ^ [k 2p k 2 + k Xll k 2 + k^kfj [Ap v ] 
+ \g»u [huad (M 2 )] - (h - k 2 f [7 log (M 2 )] 

+ i {2k 2l/ k 2fl + k lu k 2fl + kxykiv + 2k Xv k lp ) [ii og (M 2 )] 
+^(4tt)- 2 [(fci - A; 2 ) M (fci - fc 2 )„ - ^ (h - k 2 f 

]z ((h - k 2 f ■ M 2 ) - Z 2 ((h - k 2 f ■ M 2 ) 

1 



x 



- % {A 7[ )- 2 {k 1 + k 2 ) li {k 1 + k 2 ) v 



-Zo^h-hf-M 2 ) 



(31) 



(32) 



(33) 
(34) 



(35) 



In these, the functions Z (q ; M ) and Z 2 (q ; M ) are finite and can be written generically 

as 



Z k {q 2 ;M 2 ) 



/ cfe z k log 


"g 2 z(l — 2) — M 2 ' 


—M 2 


/o 



(36) 



and Dq^, A^, V^, J/ 09 and are momentum- independent divergent quantities given 
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by 



□ 



f d A k 2Ak IM k„k a kp 
J a (2tt) 4 (k 2 - M 2 ) 4 ~ 
' 4 k Akpk^ 



9a.fi 



J a (2ir 

r d 4 k 



a (2tt) 4 (k 2 - M 2 ) 3 

d 4 £; 4^ 



(2tt) 4 (jfe2 - M 2 ) 3 J a (2tt) 4 (jfc2 - M 2 ) 3 ' 

d 4 k 



d 4 k Ak^k„ 



9v* 



[XV 



(2tt) 4 {k 2 - M 2 ) 3 J a (2tt) 4 (P - M 2 ) 2 ' 
d /c 2/ci//c^, f d k fiv 



a (2tt) 4 (k 2 - M 2 ) 2 A(2vr) 4 (A; 2 -M 2 )' 



7a (2tt) 



a (2tt) 4 (P - M 2 ) 
d 4 fc 1 



4* (M 2 ) 

V4 M ') - y A J^fJW-M 2 ) 

Since the one-point functions are purely divergent, they can be expressed entirely 
of (a subset of) the above divergent quantities as: 

T s (h) = 4M { [l quad {M 2 )] + k«k{ [A/3a] | , 

T^(h) = 4{-fcf[v fr ]-^fefcnn a/)F ] 

+ l -k\kl [A^] + l^fcffcf [A Q/3 ]}. 
The two-point functions contain finite and divergent parts, and can be written as: 

T ss (h, k 2 ) = 4 j [l quad (M 2 )} + \ [AM 2 - (h - k 2 ) 2 ] \I log {M 2 )} 



-2(4vr)- 2 [4M 2 - (h - k 2 ) 2 ] 
+(k 1 + k 2 ) a (k 1 + k 2 y[A af3 ], 



-Z ((^-A; 2 ) 2 ,M 2 ) 



T FF (k u k 2 ) = A\ - [l quad (M 2 )} + - (k, - k 2 ) 2 [I log (M 2 )} 



1 



-i(Att)- 2 (k, - k 2 ) 2 

— (ki - k 2 ) a {ki - k 2 Y [A a/3 ] 
-(h + k 2 ) a (h + k 2 f [A af3 ] , 



-Zq ((h - k 2 ) , M 2 ) 
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Tl A (fa, fa) = AM(k\ - k 2 ),{[I log (M 2 )] - *(4tt)- 2 [Z ((fa - k 2 f , M 2 )] } , (46) 
-4M(h + k 2 f[A^ (47) 
-2te^ a/3 (k 2 - hfih + k 2 f [A£] , (48) 



T™(fa,k 2 ) 
T^(fa,k 2 ) 



T V J (fa, k 2 ) = i[{k x - k 2 fg, u - (h - k 2 ),(k, - k 2 ) v ] { [l log (M 2 )] 

"1 (2M 2 + (h-k 2 ^ 



-i(47r)" 



3 + (h - k 2 ) 2 



[Zq ((fa — fa) , M 2 )] 



+ V,(49) 



[(fci - fc 2 ) V - (A* - fc 2 ) /1 (A;i - fa),] { [/ W (M 2 )] 



3 

-i(4vr)- 2 



1 (2M 2 + (jfex - A; 2 



+ 



[Z ((^-fc 2 ) 2 ,M 2 )] 



3 ' (h-k 2 y 

8M 2 g^{[I log (M 2 )} - t(4n)- 2 [Z ((fa - A; 2 ) 2 , M 2 )] } + A, v . 



(50) 



T^(fa,k 2 ) 
T* T (fa,k 2 ) 



-2 { (fci - fa)^ (h + k 2 f [Afr] - (h - k 2 ) v (h + k 2 f [A^] } , (51) 



= -2ie^ al3 (ki - k 2 ) a (k x + k 2 ) 



A? 



(52) 



T^Ah,k 2 ) 
T^(h,k 2 ) 



AiMe^ih + k 2 f 



4M {g a u9/j.x - g a n9v\) (fa~fa) X 

x { [J W (M 2 )] - 2(4tt)- 2 [Zq ((fa - k 2 f , M 2 )] } , 



(53) 



(54) 



Totfav (fa, fa) 



(dafid^duX — QavQpidiiX + 9fiv9ai9^\ ~ 9Pii9ai9v\) 

x ^ [(h - k 2 f -(fa -faf (h — fa) x ] 



x <^ Ii og (m 2 ) - i(47r) 



-2 



1 2M 2 + (fci - fc 2 ) 2 rv , 2 ^ 
3 + {kl -k 2 ) 2 ^(M,(^-^ )] 



+4 - <7a^) |- [W*(M 2 )] + - [4M 2 + (fci - A; 2 ) 2 ] [J W (M 2 )] 

_I [4M 2 + (h - k 2 ) 2 ] t(4n)- 2 [Z ((fa - k 2 f , M 2 )} 
-\(fa - fa)\fa - faf [A A? ] - \(fa + fa) x (fa + faf [A A? ]| 



+9anAp v — g au Aj3^ + gp v A aix — gp^A ail . 



(55) 
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In the above, we have defined A„ v as 



4[V M „] + (ki - k 2 ) a (h - k 2 f 

x 



1 1 A 2 A 

jDa^ + 7^9auA^p + g ail Ap v — g^ u A a p — -g a pA^ u 



+ [(A* - + fca)" - (fci + fc 2 ) Q (&i - k 2 f] 



x 



1 1 1 



+(h + k 2 ) a (h + k 2 f [n a/3llu - g^A ua - g aii Ap v - 3^A a/3 ] 



(56) 



The remaining one- and two-point functions not considered above are identically zero due 
to vanishing Dirac traces. 

At this point it is important to note the generality of the method. No momentum shifts 
were done and no single divergent integral was calculated. For this reason, the results 
obtained can be used with any preferred regularization method. Of course, different reg- 
ularization schemes can lead to ambiguities and violation of symmetry relations between 
Green's functions, as we shall discuss in the next section. 



IV. AMBIGUITIES AND THE CONSISTENCY CONSTRAINTS 

In the preceding section we have completed the evaluation of all fermionic one- and two- 
point amplitudes. Before any further steps in evaluating these amplitudes, it is important 
to notice that there are ambiguities associated with the choices for the momentum routing 
in the internal lines of loops. In the calculations performed in the previous section, we have 
not made any momentum shifts in intermediary steps and have left the labels of the internal 
momenta k\ and k 2 completely arbitrary and unspecified. However, due to momentum 
conservation a physical amplitude can depend only on the difference q = k± — k 2 , and 
no dependence on the sum Q = k\ + k 2 can be present in the amplitude. There can be 
no dependence on Q because this combination is ambiguous: two different choices of k\ 
and k 2 can give the same q, but they will always give different Q's. In view of this, one 
could imagine that it is impossible to make any predictions because the amplitudes contain 
dependencies on arbitrary contributions which are functions of Q. In general, depending on 
the regularization scheme used, this is actually true. In the literature associated with the 
NJL model, it is usual to use a particular routing, the symmetric combination k\ — q/2 and 
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&2 = —9/2, meaning that Q = and the amplitudes are automatically free of ambiguities. 
However, restriction to one particular choice, besides breaking homogeneity of space-time, 
also leads to difficulties with amplitudes containing more than two propagators where there 
are in principle more than two arbitrary momenta and choices like the above would not be 
allowed because in general would lead to violations of momentum conservation. 

Our aim now is to show that what we have done so far does not lead to violations of the 
relations among Green's functions because of arbitrary choices of the momenta k\ and &2- 
As shown in Section [Til such relations are very general, but due to the ultraviolet divergences 
they can very easily be violated when not being careful with the explicit evaluation of the 
integrals within a particular regularization scheme. 

It is a simple task to identify the ambiguous terms, proportional to Q. In the one-point 
functions they are given by (note that for these Q = k\): 




(57) 



(58) 



In the two-point functions, the ambiguous terms are given by: 




(59) 
(60) 
(61) 
(62) 







(63) 





+Q a Q po^i, — g^pAvu — g ail lS.$ v — 3^A a/ 3] , 



(64) 
(65) 




(66) 
(67) 
(68) 
(69) 
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+9P» { T *») ambi ~ 9PH ( T av) ambi (70) 

Note that all ambiguous terms appear as coefficients of the divergent quantities 
V^, and A a p. These quantities also play an important role in the analysis of 

the algebraic relations involving different Green's functions, the consistency constrains ob- 
tained in Section [Til Those constraints were obtained by making the only assumption of the 
validity of the linearity of the integration operation. Although this assumption seems very 
reasonable, the mathematical indefiniteness due to the divergences turns the preservation of 
such relations a non trivial supposition. 

Let us now show that despite the amplitudes themselves contain ambiguous terms, the 
general relations involving different Green's functions obtained in Section [TT] are not violated 
by the manipulations done so far. We start with the relations in Eqs. ([7]) and ([8]) involving 
the Green's functions TYj ik\, k 2 ) and Tj 7 ik\). Contracting TYj (ki, k 2 ) given in Eq. (]49l) 
with the external momentum q leads to 

er^{k 1 M) = <tA lw . (71) 

From Eq. (jHJ), one sees that one must identify on the r.h.s. of Eq. ( jTTi) the difference between 
two vector one-point functions corresponding to the internal propagator carrying arbitrary 
momenta k% and k 2 . For this purpose we take Eq. (1561) for A^ u and contract it with 

q»A ia , = T?(k 2 )-T?(k 1 ), (72) 

and also with q v 

q"A tlu = T^(k 2 )-T^(k 1 ), (73) 

where we used Eq. f H3|) for identifying T ; ] / (fc 1 ). These results imply that the algebraic 
relations of Eqs. ([7]) and ([5]) involving {k\) and T^f (ki, k 2 ) obtained by the formal 
manipulations in Section [TTJ are preserved by the explicit and independent evaluations of 
(ki, k 2 ) and (ki) in Section [TTTI Therefore, the relations are preserved in spite of 
the fact that both amplitudes T„ (ki) and T^f (ki, k 2 ) have ambiguous pieces, which are 
proportional to the divergent quantities V^, n a/3fll/ and A Q/3 . 
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Next we consider the relation given in Eq. (jTTI) involving the amplitudes T^ s and T s . 
Using the explicit expression for T^ s (ki, k 2 ) given in Eq. (jlTj) , one obtains 

q^ s (k h k 2 ) = -AM (k»k\ + Kk\ - k»k\ - [AJ. (74) 

Due to the obvious symmetry under interchange of the Lorentz indexes £ and \i in A^, and 
by comparing with the explicit result for T s in Eq. (jUD , it is very simple to show that the 
identity of Eq. (1111) is also preserved. 

Now we consider the identity given in Eq. ( ITUl) . involving the amplitudes T^ 4 , T AP and 
Tp. Contracting T^f given in Eq. fl50l with g M , we get 

<f T// (fei, fc 2 ) = 2M {-4Mg, [j Zofl (M 2 ) - i(Air)- 2 Z (g 2 ; M 2 )] } + gM^. (75) 

It is easy to verify that using Eqs. and (JI2D, the identity in Eq. (jTUJ) is preserved. 

We proceed examining the relation given in Eq. (j!2p using Eq. (1461) for the amplitude 
T AP (k±, k 2 ). Adding and subtracting scalar one-point functions carrying momenta k\ and 
k 2 , one can easily check that this expression for T pp (fc 1? k 2 ) naturally leads to preservation of 
the relation in Eq. (Jl2i) . Also, the relations involving the AV two-point function in Eqs. (JT3l) 
and (|14l) are immediate: the Green's functions T PV and and T A on the r.h.s. of Eqs. (TO?]) 
and (1T41) are identically zero due to properties the Dirac traces, and on the l.h.s. one also 
obtains a zero because the contractions of the explicit expression for T^f (kx,k 2 ) in Eq. flU} 
with {ki — k 2 Y or {ki — k 2 y L are zero due to the antisymmetry of e^ v \^. 

Next, we turn our attention to the tensorial amplitudes. We start with the relations of 
Eqs. ( 1201) and ( 1251) . Initially we note that the amplitude T^J V (ki, k 2 ) given in Eq. ( 1541) can 
be reorganized as 

TZ **) = ^ (A:i ' ^ " T - ^' ^ 

" 2^ ^ PC **) - (k 1} k 2 )] , (76) 

where we have used Eqs. (j4"9l and (1501) . Written in this form, it is now trivial to see that 
Eqs. (1201) and (1251) are satisfied. Next we consider the relation in Eq. (fT8l) . The TT two-point 
function, given by Eq. (1551) . after some algebraic effort can be put in the form 

Tafav ^2) = 9 af iTp,f (ki, k 2 ) — g au Tpf (hi, k 2 ) 
+g^ A (h, k 2 ) - g^Ti A (A*, fe) 
+ (9 a »9pu- g a u9^)T ss (k h k 2 ) , (77) 
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where the Eqs.( l50p and (1441) have been used. It is now evident that the contraction with the 
external momentum leads to Eq. ([TBI . Now we consider the constraint given in Eq. fl24l) . 
First we note that Eq. (|5"Tj) for the ST amplitude, by using Eqs. (14"71) and (TTo^) . can be 
rewritten as 

T S J fc 2 ) = fcrr (fa, fa) - q v T s / (fa, fa)] , 

= -^(T^fi^h). (78) 

Therefore, contraction with the external momentum 

q"T T J (*!, = ^ **) - [T s (fa) - T s (fa)] } , 

= ~e^(T AV ) X Hh,k 2 ), (79) 

leads immediately to Eq. (1241) . Now, noting that by Eqs. (|4"8|) and (|52|) the Ay Green's 
function is identical to the TP function, we have that 

q^f (fa, fa) = -2MT^ P (fa, fa) , (80) 

as it should to satisfy Eq. ( I2TI) . Also, comparing Eqs. ( 1471) and ( 1531) . one sees that 

^ fe) = ~<W (T SV f fc) , (81) 

and then 

(fa, fa) = 2MT£ P (fa, fa) , (82) 

q»T™ (fa, fa) = -2MT™ (fa, fa) = -z^g A {T sv f (fa, fa) . (83) 

Therefore, the remaining relations, given by Eqs. ( 1221) and ( 1231) . are also satisfied. 

This completes the verification of the consistency of the manipulations performed at the 
one loop level. It is important to emphasize that all the relations among Green's functions 
which can be stated at the level of integrands are preserved, in spite of the presence of am- 
biguous terms. The arbitrariness concerning the choice of the regularization is also preserved 
since only very general mathematical properties have been assumed. 

The crucial point now is that the preservation of the relations among Green's functions 
is not the only requirement one should ask for a consistent regularization scheme since, 
as we have seen, the amplitudes themselves contain ambiguities that are functions of the 
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ambiguous combination of momenta Q = k\ + k%. Therefore, one must have a scheme that 
ensures elimination of the ambiguous terms from the amplitudes themselves, not only in 
the relations involving two or more of them. In the next section we discuss the constraints 
imposed by symmetry relations (Furry's theorem and Ward identities) on special amplitudes 
and they will provide guidance for dealing with the ambiguous terms. 

V. AMBIGUITIES AND PRESERVATION OF SYMMETRIES 

In the process of constructing a consistent interpretation for the divergent one-loop am- 
plitudes the preservation of symmetries plays a central role. In principle, there is no a priori 
reason for expecting that space-time symmetries will be automatically manifest in divergent 
amplitudes. However, it seems nevertheless reasonable to expect that one should be able to 
identify general properties that the divergent quantities must satisfy in order to guarantee 
the preservation of such fundamental symmetries. In this sense, the regularization method 
itself is not the most important ingredient, what really matters are the requirements that 
quantities like V^, n a/3fll/ and A a p must obey to preserve the symmetry relations. Having 
this in mind let us now consider the symmetry properties pertinent to the one- and two- 
point amplitudes we are discussing. We shall refer to the Ward identities and other general 
constraints imposed by Furry's theorem on these amplitudes. 

We start considering the simplest amplitude that carries one Lorentz vector index, the 
amplitude (ki). On general symmetry grounds, Furry's theorem states that this ampli- 
tude must be zero. So, from Eq. f|43l) . Furry's theorem requires that 



There are two different ways to satisfy this requirement. The first one is the choice k\ = 0. 
But, is it always possible to make this choice? Thinking on in isolation, the answer to 
this question is affirmative, since k\ is arbitrary. However, is not the only amplitude 
in the theory and so one must ask the question if this choice is not invalidating other 
symmetry relations. For example, the relation given in Eq. (JH1) relates this amplitude to 
(ki — k'lfT^ (ki, k2)i which we repeat here for clarity, is given by 




•v 






(85) 
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Vector current conservation demands that (k\ — k^fT^ (k\, k 2 ) = 0. Therefore, the differ- 
ence of the two one-point functions on the r.h.s. of Eq. (185]) having dependencies on k\ and 
k 2 needs to be zero. Obviously, the simultaneous choice ki = and k 2 = 0, which would 
satisfy both requirements, cannot be made because this would imply q = k\ — k 2 = al- 
ways. Therefore, we need another way to satisfy Eq. (I84p . Since the requirement of Eq. (I84p 
involves the divergent quantities O a pp V , V ^ and A Ml/ , one could ask for a regularization 
scheme that leads to 

□3* = Kl 9 = = °> (86) 

where the superscript reg means that the integrals defining these quantities are regularized. 

The same conclusion is reached considering the explicit expression for TYj {ki,k 2 ) given 
in Eq. ( 141)]) . Contracting it with the external momentum g M = (k\ — k 2 )^, one obtains: 

<?r™ (h, h) = 4 L a [v w ] + - k a 2 k^k p 2 ) 1 - [n aPpu ] 

-{klk{ - k\k p 2 ) X - [A pv ] - {k lv k a x k{ - k 2v k*k p 2 ) 2 - [A a/3 ] J . (87) 

Since a conserved vector current should not be obtained by convenient choices of the arbitrary 
momenta k\ and k 2 , the conditions of Eq. (}86j) seem therefore also necessary here. 

For the same reason that (fci) must vanish, other vector two-point functions need 
vanish identically. These are T^ s (ki, k 2 ) and T^f (fci, k 2 ) which, from Eqs. (l4"Tj) and (l4"8l . 
imply in 

Q« [A^] = 0, (88) 
e^p^Qt [A« Q ] = 0. (89) 

In principle, for these two specific amplitudes both options, of choosing k\ and k 2 in a conve- 
nient way or constructing A r ^ = 0, are possible. For example, considering the contractions 
of these amplitudes with external momenta we obtain 

q^ s (k u k 2 ) = —AMq^Q 13 [A^] , (90) 
<?T™ (k u k 2 ) = -2e, va pq^Q i [A« a ] , (91) 
fT™ (h, k 2 ) = -te^fqfQt [A* a ] . (92) 

A conserved vector current for T^ s (ki,k 2 ) can be obtained with the choice ki = —k 2 , or 
by taking A r ^f = 0. However, both contractions involving T^f (k±, k 2 ) vanish identically 
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independently of the two possible choices, just because the antisymmetric e^ va p is contracted 
with a symmetric object. The vector current must be conserved, but the axial- vector current 
must not. So there is only one consistent value for T A „ (ki, k 2 ): the identically zero value. 
Otherwise a symmetry relation is broken. We can add to this argumentation another very 
general aspect that forces us to obtain a zero value for T A ^ (ki, k 2 ) (and T^ s ): unitarity. 
If the amplitude does not vanish then it needs to develop an imaginary part at q 2 = 4M 2 
to be consistent with unitarity (Cutkosky's rules). Clearly, from Eqs. fj47|) and fj48|) for 
T^ s (ki, k 2 ) and T A ^f (ki, k 2 ), respectively, this cannot happen. 

Next we consider T A (ki,k 2 ). Using its explicit expression given in Eq. ( l50i) . one can 
show that 



However, the proportionality between the axial-vector and the pseudoscalar current states 
that q»T AA (k x ,k 2 ) = -2MiT? A (k u k 2 ). Therefore, one arrives at the same conclusion as 
for the amplitude T^f (k%, k 2 ), that the relations given in Eq. (1861) must be satisfied, since 
T^ A (k\, k 2 ) is free from ambiguities - see Eq. ( T4*6l) . Also, the same conclusion is obtained if 
one rewrites q^T AA in terms of the amplitudes T^f (k\), as in Eq. ( TlOl) . 

Considering all amplitudes and their symmetry relations, the same conditions will emerge: 
there is no consistent interpretation for the one-loop divergent amplitudes if the conditions 
given in Eq. (1861) . that we call consistency relations (CR), are not fulfilled. In principle, one 
could argue that the imposition of the CR's represents an arbitrary choice which is at the 
same level of the choice of a specific regulating distribution to be used in the integrands 
of the divergent integrals. However, this is not true. One should consider the CR's as a 
fundamental requirement to be imposed on the one-loop divergent amplitudes in order to 
materialize the fundamental space-time symmetries. Any calculation that violates the CR's 
has the potential of predicting unphysical results, since it leads to the destruction of the 
foundations of the theory which have generated the amplitudes themselves. Therefore, the 
CR's do not represent arbitrary choices, because there is no option out of these properties 
capable to allow a consistent interpretation of the calculations. Note that the CR's not only 



q"T AA (h, k 2 ) = 4 j-g" [V au ] + {k^k p x - k^k^ - [U aPpv ] 

~\~ {k\k^ k 2 k!p) — + (^ki u k^ k^ k 2u k 2 k 2 
-2MiT? A {k u k 2 ). 



(93) 
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remove all the ambiguous terms, which are always symmetry-violating, but also remove all 
the symmetry violating terms, which are not always ambiguous. 

To finalize this section, we summarize results by defining what we denominate the "con- 
sistently regularized amplitudes " (CRA) denoted by T s , T^ s , T ss , • ■ ■ . These are respec- 
tively the amplitudes T s ,T^ S ,T SS , ■ ■ ■ obtained previously, with the terms containing the 
pieces proportional to the quantities Da^, V Mi ,, and removed, as demanded by the 
CR's given in Eq. fl86l) . Explicitly, they are given by: 

I) One point functions: 



T s = AM [l quad (M 2 )] , 



V = 0. 



(94) 
(95) 



II) Two point functions: 



(96) 



r SS (g) = 4<( [I quad (Af 2 )] + i (4M 2 - g 2 ) [l log (M 2 )] 



-- (AM 2 -q*)(^- 
2 V H ' \ 16tt 2 



[Z (q 2 ,M 2 )] 



(97) 



T pp (q) = A{-[l quad {M 2 )] + \q 2 [l log {M 2 )] 



I 2( 1 

2 q V16tt 2 



[Z (q 2 ,M 2 )] 



(98) 



T PA (q)=AMqA[l log (M 2 )]- 



16tt 2 



[Zo(q 2 ,M 2 )] 



(99) 



[llog(M' 



167T 2 



(100) 



r* A (?) 



g (?^M^ - ?A»9") 



X 



[«M 2 )] 



167T 1 



-8MV J [//o,(M 2 )] 



I + (^1±«!) [W ,M*)] 
z 



16tt 5 



[^o(g 2 ,M 2 )] 



101) 
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The tensor amplitudes can be written as 



T?(q) = T? T (q) = T a fAd) = 0, (102) 
(q) = AM (g au g, x - g a ,g uX ) q x { [I log (M 2 )] - [Z (q 2 , M 2 )] J , (103) 



Taflfiv (q) = (9an9/3tgv\ - 9au9t3(9fi\ + 9Hv9a£,9ii\ ~ 9^9 0^9 , 



4 (A £A - aV) {l, 3 (M 2 ) - (J^j [i + [Z„ (A/ 2 , q 2 )] 

+4 (vsj, - g w 9*) |- [l,uad(M 2 )] + i (4A/ 2 + j 2 ) [/ IO! ,(A/' 2 )] 

^(4M 2 + 9 2 )( T ^)[Z„(, 2 ,M 2 )]}. (104) 



VI. PHENOMENOLOGY AND NUMERICAL RESULTS 

The amplitudes obtained in the preceding section are free of ambiguities and are sym- 
metry preserving. In this paper we have focused on one and two-point functions, where 
resides the highest degree of divergence, however similar results can be obtained for three 
and four-point functions used in the model to describe meson decays and meson-meson in- 
teractions. All these mathematical structures also appear in fundamental theories. Here, 
as a consequence of the adopted strategy to handle the divergent structures, it must be 
noted that such structures are treated in a very closely related way as in renormalization 
procedures of renormalizable theories. In the obtained expressions for the calculated ampli- 
tudes only two divergent objects have survived after the adoption of CR's, namely Ii og (M 2 ) 
and I qua d(M 2 ). The next step, if the amplitudes were to be considered in the context of 
fundamental theories, is the elimination of such objects through the reparametrization of 
the theory in the renormalization of physical parameters. Since Ii og (M 2 ) and I qua d(M 2 ) are 
completely absorbed in this process, the regularization eventually used plays no relevant 
role, due to the fact that the renormalized expressions are independent of particular aspects 
of the chosen regularization. The theory is predictive, given the fact that the results for the 
amplitudes associated to physical processes do not depend on the choices involved in the 
intermediary steps. 

In the case of NJL model, considering the nonrenormalizable character, the remaining 
undefined objects need to be specified in some way by using physical parameters chosen 
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as inputs of the model. The traditional way to make predictions in the context of NJL 
involves regularization such as the introduction of a 3D- or 4D-cutoff A in the involved 
Feynman integrals. The expressions for the amplitudes are, in this way, dependent on 
the parameters of the chosen regularization distribution as well as, if the CR's are not 
satisfied, on the chosen routing for internal lines momenta. In this scenario the model 
contains at least two parameters (in the chiral limit m = 0); the coupling strength Gs 
and a regularization parameter (the cutoff A). The constituent quark mass M is not an 
input parameter as it is given by the solution of the gap equation. The parameters Gs 
and A need to be fixed through two experimental information. The quark condensate (tpip) 
and the pion decay constant f n are usually used for this purpose. The first is related to a 
quadratically divergent amplitude while the second is related to a logarithmically divergent 
one. The values for Gs and A which give the best adjustments to the experimental values 
of (ipip) and f n will depend on the specific form of the regularization distribution. As a 
consequence all the physical amplitudes describing processes pertinent to the model are 
affected by the choice of the regularization. Due to this reason the chosen regularization 
must be part of the model. The predictive power of the original quantum field model is 
affected since the predictions are dependent on a choice which characterizes an ambiguity. 
This is not the desirable situation. We wish any model prediction becomes unique in a 
similar way as it happens in renormalizable theories. The regularization must become just 
a convenient choice in the intermediary steps. In order to show that it is possible to achieve 
the desirable situation, referred above, we first note that the adoption of the CRA's, listed 
in the preceding Section, implies that we can only adopt regularizations that fulfill the CR's. 
This does not represent a choice as, following our analysis, regularizations which break the 
CR's will lead to ambiguities as well as to symmetry violations and this is unacceptable. If 
we want a predictive model, only regularizations preserving the CR's make sense . After 
these important remarks let us now show how the arbitrariness associated to the choice of 
regularization can be completely removed. 

Having this in mind we start by showing how to relate the remaining objects Ii og (M 2 ) 
and I qua d{M 2 ) to physical observables chosen as input of the model. First we point out that 
the quark condensate (ipip) is related to the T s one-point function as 

(faj>) = -N C T S . (105) 
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Substituting the result fl94l) we get 

= —AN C M [il quad (M 2 )] . (106) 

For the second we note that, in the context of NJL model, mesons are relativistic quark- 
antiquark bound states. In the random-phase approximation (RPA), the meson propagators 
can be written as (see for example Ref. |5j) 

D » W - ; _ 2 a% M m ■ m 

where Um is the polarization function defined by 

/d A k 
— -^Tr {T M S (k + ki) T M S (k + k 2 )}, (108) 
(z7r) 

with S being the quark propagator defined previously. Tm stands for the flavor and Dirac 
matrices giving the quantum numbers of the meson A4. For example, for the neutral pion, 
T M = r 3 7 5 , for the scalar-isoscalar meson, Y M = 1. In writing the equations above, we 
assumed the most general labels for the momenta k\ and k 2 running in the internal lines of 
the loop integral. The physical momentum q is defined as the difference k\ — k 2 as imposed 
by energy-momentum conservation at each vertex. 

The pole of the propagator in Eq. (11071) . calculated at q 2 = m 2 M) gives the mass of the 
respective meson. The condition for the pion mass is given by 

1 - 2G s IL K {ml) = , (109) 

where 

U n (q 2 ) = -iN c N f [T pp (q)] . (110) 

Substituting now the explicit form of T pp (q), derived previously and given in Eq. (!98|) . we 
obtain 

(q 2 ) = 4iN c N f [I guad (M 2 )} 

- 2 iNc N f q 2 { [l log (M 2 )] - (J^j [Z {q\ M 2 )] } . (Ill) 

Using the Eq. d3J) in order to eliminate I qua d (M 2 ) and evaluating Eq. (11 111) at q 2 = m 2 we 
get the following expression for the pion mass 

m 2 = ^ I (112) 

^ AN c N f MG s [I lo9 (M 2 )] + (^) [Z (ml, M 2 )]} ' 1 ) 
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As seen, in the chiral limit (mo = 0) the pion becomes massless (m n = 0), in agreement with 
Goldstone's theorem. 

The pion phenomenology is also characterized by the decay constant /„.. Experimentally 
it is related to the weak decay it^ —>■ ^ + and is calculated from the vacuum to one-pion 
axial- vector current matrix element 

(01^)7^/2^)1^' (g)) = i Uq^e^*, (113) 

where \^{q)) is a pion state with four- momentum q. At one- loop order, one can express 
this matrix element in terms of the T^ p two-point function as 

i/ w <7A = -N c g nqq 6 l3 [Tf p (q)] , (114) 
where g nqq is the pion-to-quark-quark coupling strength, related to the residue of Eq. fllQTf) 

dU PP (g 2 )V X 



as 



2 

9nqq ~ I oo 



(115) 



dq 

Using Eqs. (}99l) and fill II) . we can write 

U = -AN c9nqq M |z [l log (M 2 )] + (jL^j [Z (ml M 2 )] } , (116) 
9% = —2N c Nf [l log (M 2 )] + (j^j [Z K,M 2 )] } 

- 2N c N f ml (jL^j [Y, (ml, M 2 )] (117) 
where Y\ (q 2 , M 2 ) is the k = 1 element of the set 

Y k (q 2 ;m 2 ) = fdz - Zk{1 - Z) - . (118) 

Jo q z (i — z ) — m 

In a completely similar way, for the scalar meson (a) we have 

m 2 =4M 2 - — - (119) 

MiG s N c N f {iI log (M 2 ) + (j^)[Z ( m l,M 2 )}Y 1 ' 



g~ 2 q = 2iN c N f (ml - AM 2 ) \ [Y, (m 2 , M 2 )] 

- 2iN c N f |/ W (M 2 ) - (r^j [Z Q K , M 2 )] J . (120) 

It is easy to see from Eqs. f 1 11 6 j) and (11171) that the observable f n can be related to the 
undefined quantity Ii og . In order to make the aspects we want to emphasize clear, we initially 
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consider this relation in the chirally symmetric case. Then we get 



a >°> ( M2 ) = -2mw (121) 

since Zo(m% = 0) = 0. Through the Eqs. (I106P and (j!2ip we have stated relations between 
two observables (inputs), the quark condensate and the pion decay constant, with two unde- 
fined objects coming from loop calculations. These objects are functions of the constituent 
quark mass which must be determined in some stage. Therefore the Eqs. (I106p and f ll21j) 
in fact represent the relation of two quantities to two functions. If we want to know such 
functions we have to integrate Ii og (M 2 ) and I qua d(M 2 ) which means to adopt an explicit form 
of regularization distribution. This process introduces at least one regularization parameter 
A as it is well known. Different regularizations will generally lead to different values of A as 
it is usual in the context of NJL model with regularizations. In order to avoid this situation 
we will proceed in a different way. 

First we note that those two functions are not independent. It is possible to show that 
they are related by 

[tI quad (M 2 )]=tI log (M 2 ). (122) 



DM 2 

On the other hand, we can also state that Ji og (M 2 ) possesses the following property 

^[ i7 '°*( M2 )] = l» (123) 
In order to satisfy these two conditions it is necessary to get the following general forms 

tI log {M 2 ) =^lnM 2 + C 1 , (124) 
i huad (M 2 ) = -^— 2 M 2 [In M 2 - 1 + Ci] + C 2 , (125) 

where C\ and Ci are indeterminate constants - C\ is dimensionless and C% has dimension 
of (mass) 2 . In the context of regularizations, C% and C2 are related to the regularization 
parameter A. Eliminating the constant C\ we see that 

i huad (M 2 ) = + M 2 [i J log (M 2 )] + C 2 . (126) 

Replacing I qua d and Ii og in term of (ipip) and /„-, Eqs. (11211) and (11061) . we get 

+ -Jjl CAM- ±— -L = 0. (127) 



16tt 2 \2N c N f J AN, 
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There are two important aspects involved in the above equation. First, it ensures that it is 
crucial to obey the properties (11221) and (11231) when the functions Ii og (M 2 ) and I qua d{M 2 ) 
are made explicit. These two properties work as additional constraints to be required of a 
regularization distribution if one wants to get consistency in perturbative calculations. The 
violation of these properties will result in breaking the scale properties of physical ampli- 
tudes, and possesses the same status of symmetry violations 25[]. The second aspect refers 
to the dependence of the physical parameters on the choice of the specific regularization. 
The above expression states that even if a regularization obeys the CR's and simultaneously 
the conditions H122I) and (I123p it remains a freedom to distinguish it from other regular- 
izations belonging to the class of consistent regularizations, which is the value of Ci- Two 
consistent regularizations can differ only by the value of C<i- As a consequence, the physical 
implications of the model seem to be definitely regularization dependent. The introduction 
of experimental values for the inputs (ipip) an d fn, it makes necessary to specify C2 to get 
M and then through the gap equation, to get the value for Gs- There is nothing more to 
be imposed, based on consistency reasons, to remove this arbitrariness. Apparently all our 
efforts cannot avoid the dependence of the results on the choice of the regularization, even 
if we have drastically restricted the regularizations which can be used in the calculations. 

At this point it seems there is nothing else to do except to choose convenient values for 
C*2- The convenience of such choices is related to the fact that we need to choose an adequate 
value for C2 in order to get a good value for M and after this to find a good value for Gs- It 
is precisely in this process that emerges the most surprising aspect of our formulation. If we 
recognize that only positive values of M make sense in the equation, the nonlinear character 
of the equation in M produces a critical condition to possible values of G^- It turns out that 
one finds solutions with M > only for C2 > C cr it, such that 

C2 < C cr u — > no solutions 

C2 = C cr it — ► one solution (128) 
C2 > Ccrit two solutions. 
Therefore, it seems obvious that there is only one value for C2 which is reasonable, the 
Ccrit, due to the fact that only this value allows us a consistent physical interpretation of 
the model predictions. Assuming this attitude, the NJL model becomes predictive since 
the remaining arbitrariness is fixed through the existence of a critical condition. There are 
considerable differences between our formulation and the traditional ones. The first and 
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immediate refers to the determination of the constituent quark mass M. Such value is fixed 
by the critical point in the diagram M x C%. The value of the mass, therefore, depends only 
on the experimental values for (ipip) and f w which are the chosen inputs of the model. The 
gap equation will be used in the determination of the coupling G,s compatible with the mass 
fixed by the critical condition. At this point it is crucial to ask: are the values for M and 
Gs, emerging from this critical condition, reasonable? 

In Fig. [1] we plot all possible physical solutions (M > 0) of Eq. (I127j) as function of the 
arbitrary constant C^- Using as input (^V) = (—250.0 MeV) 3 and f n = 93.0 MeV it results 
that for values of C2 < C cr u, Eq. (11 2 7ft there is no physical solution, at C2 = C cr i t there 
is only one solution, and for C2 > C cr u there are two possible solutions. In particular, at 
the critical point we obtain C crit ~ 24.82 MeV 2 , M ~ 468.4 MeV and G s ^ 7.5 GeV~ 2 . 
Therefore the values for M and Gs are in good agreement with the ones found in literature 
of this issue. 

800 

> 600 
5 
5 

400 



"24 27 30 33 36 

C 2 (MeV) 2 

FIG. 1: Solutions of Eq, (|127p as function of the constant C<i- 

Another difference between this formulation and the traditional ones refers to the meson 
phenomenology. When we adopt the present formulation, in the presence of chiral symmetry 
breaking parameter m , previously fixed, since it is an input parameter, the meson masses 
and meson- quark- quark couplings, as well as other physical aspects, emerge as genuine 
predictions, including those of pion. In order to see this aspect clearly let us consider the 
case with rrtQ 7^ 0. As it can be seen in Eq. f)112p . the pion mass is nonzero, reflecting the 
fact that the original Lagrangian is not chirally symmetric with itlq 7^ 0. The introduction 
of a current quark mass modifies the results for other quantities as well. For example, the 
expression for f n now contains a finite part, and because of this, the expression of Ii og (M 2 ) 
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in terms of f 2 changes to 

[Z (ml, M 2 )] 



* [lio 9 (M 2 )] 



16tt 2 



ft L / N c NfmlM 2 \Yi (ml, M 2 )] 
AN c NfM 2 { + V 2tt 2 / 2 



(129) 



As a result, we can write for the pion mass the following expression 



2 _m M N c N f mlM 2 [Y, (ml, M 2 )] , 

Also, for the effective pion-coupling constant we obtain 
9^ = -^^ [Y,{ml,M 2 )] 



ft L, L N c N f m 2 w M 2 [Yi (ml, M 2 )} \ 
+ 2M 2 { 1 + f 2^71 J' (3) 

As before, eliminating I quac i(M 2 ) and Ii og (M 2 ) in favor of (^pip) and f n , we obtain a nonlinear 
equation for M that now involves the pion mass m n , 

£ [1 + * WlJO] _M + ft M 

This equation is much more complicated to solve the one in the case of exact symmetry. 
This equation and the expression for m n , however, can be simplified using the following 
approximations for Z and Y\ 

Z (q 2 - m 2 ) = - + • • • , (133) 

KH ' ; 6m 2 60m 4 ' V ; 



and 



M^) = ~-^ t+ -- (134) 



Using these approximations, we obtain 

, m„M / / N c N f mj \" 



A ^ + ^L( 1 + 1 / 1 + Mf£\ _ <E _ c 2M -M = 0. (136) 



16tt 2 AN c N f V V 127r ^- / 967r2 4iVc 
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Inserting the "experimental" value for mo (as well as for and f n ) we search for the value 
of Gi which corresponds to only one positive value for M and thus determine the values for 
Gs, tn-Ki and so on. For this propose we take mo = 5.2 MeV, obtaining m n ~ 135.3 MeV, 
M ~ 468.4 MeV, G s ^ 7.5 GeV" 2 and C cnt ~ 24.82 MeV 2 . Finally, Eq.flEDJ furnish 
g nqq ~ 4.97. These predictions are in good agreement with experimental data and those 
used in the literature. 

In the scalar channel, we have the following expressions for the a meson mass and the 
aqq coupling constant 



m„ 



4M 2 - 



MGf 



N c N f 
4tt 2 



Z (ml,M 2 ) + 



ml 



QM 2 




N c N f ml 
12vr 2 / 2 

Z (m 2 ,M 2 ) + 



-i 



(137) 



6M 2 



+ 



^(m 2 -4M 2 ) [^K,! 2 )] 



2M 2 



1 + 




N c N f ml 

127T 2 / 2 



(138) 



Numerically we have m a = 938 MeV and g, 



aqq 



2.29. 



The fact of the arbitrary character of Ci being removed, owing to the existence of a 
critical condition, it is the most important result point in the analysis made in this Sec- 
tion. This means that the phenomenology becomes completely independent of the specific 
regularization scheme employed if such scheme is consistent with the scale invariance. 

The model, within the scope of this prescription, becomes predictive since the role played 
by a regularization has completely disappeared. In this sense, in spite of being a nonrenor- 
malizable model, the predictions are made in the same spirit as in renormalized models, 
at the considered level of approximation. This is, undoubtedly, a very important improve- 
ment in the quality of this type of perturbative calculations. However, we must be aware 
of the fact that this does not represent the solution to all the problems involved. Since 
the amplitudes have acquired structures which are very similar to those belonging to the 
renormalized theories, what remains to be considered are the phenomenological implications 
for the model predictions of the so called Landau vacuum instabilities or ghost poles in bo- 
son propagators. This aspect of the perturbative calculations appeared in hadron physics 
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in connection with the Skyrme model 2y, |27|, |28|, |29|]. More specifically, Landau vacuum 



instabilities were found in chiral-quark models when soliton solutions were searched for in 
the renormalized sigma model. (For an update of this problem please see the recent paper of 
Arriola-Broniowski-Golli [2^]). However, Landau instabilities seem to be present in almost 



all renormalizable theories where fermions are coupled with boson fields 



. Their occur- 



rence can be indicated by the presence of tachyon poles in boson propagators corrected by 



one-loop fermionic contributions 



3jj. In asymptotically free theories where the bosons are 



self-interacting fields, contributions coming from bosonic one-loop diagrams may eliminate 
;he problem 3J|, otherwise the general rule seems to be the existence of vacuum instabilities 



281 ] . In the QED, where the vector gauge field is not a self-interacting field, tachyonic poles 



occur in the photon propagator such that the vacuum is unstable at the one-loop level. In 
the linear sigma model the meson propagators are equally contaminated by tachyonic poles. 
In the QED, the presence of such type of undesirable poles does not play a physically rel- 
evant role because the scale of the fluctuations at which the instabilities occur is me* (m 
is the electron mass and a is the fine structure constant) 28]. Therefore, this happens in 
a region which is certainly beyond the expected validity of the theory. On the other hand, 
in hadronic phenomenological models, where a fermionic field is coupled with a mesonic 
one, like in the Yukawa model or in the chiral a model, such scale of flutuations changes 
drastically (around lGeV) due to the nucleon or quark mass and due to the value of the 
constant coupling involved. This means that the ghost poles may have relevant influence in 
phenomenological implications of the model. 

In NJL model we have only fermions in the Lagrangian but the meson states are inter- 
preted as a quark-antiquark bound states. The intermediate amplitudes are fermionic loops 
such that, in the RPA approximation, the meson mass is identified as the pole of the Eq. 
(107). This is precisely the structure of the renormalized meson propagator in the linear 
a model. This means that if we look carefully to the condition stating the pion mass, for 
example, we will find that in addition to the pion pole, there is a tachyonic pole in the 
Euclidian region with a negative residue. This implies that the corresponding dispersion 
relation will be verified only if this pole is included. This situation is not commonly consid- 
ered within the context of NJL model due to the use of cut-offs which change the behavior 
of the fermionic one-loop contributions such that the problem is automatically eradicated 

y. 
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In the procedure adopted in the present work the finite parts of the Green's functions 
are not modified, putting the physical amplitudes at the same level as those belonging to 
renormalizable theories. Because of this, the questions related to the ghost poles or vacuum 
instabilities may become relevant. The presence of ghost poles in meson propagators or 
Landau vacuum instabilities in the NJL model is expected due to its equivalence to the 
linear a model. Although the Green's functions are, strictly speaking, not the same ones, 
due to the definition of the renormalization parameters, the S'-matrices of both theories are 
identical 

On general grounds it is not completely clarified if ghost poles are real ingredients of QFT 
or if they are a product of a particular kind of perturbative solution (one-loop approximation) 



28j | . Due to the fact that if they are real ingredients of QFT fundamental axioms are violated 



since only real poles are expected to exist, the relevant question seems to investigate if 
this type of instabilities survives to higher-orders calculations. The acceptance that these 
undesirable poles are unavoidable and, therefore, real aspects of QFT, constitutes a very 
frustrating fact just because we are accepting that the solutions we can obtain do not 
obey the fundamental axioms of the theoretical apparatus we have constructed. The most 
reasonable expectation is that the instabilities will disappear when contributions of higher- 
order are computed. If this is the case, there is nothing else to do but concluding that 
the one-loop approximations are not adequate to investigate phenomenological implications 
of a theory or model. If, however, the instabilities constitute an unavoidable aspect of 
certain classes of QFT, we must make efforts to get theories free from those problems by 
construction, like in the case of anomalies. Only additional investigations will clarify these 
doubts. We are, however, convinced that the questions related to the regularizations in QFT 
are of a different nature from those related to the ghost poles in perturbative corrections of 
boson propagators. 



VII. CONCLUSIONS AND PERSPECTIVES 



We considered in detail questions relative to the predictive power of the NJL model. 
Given its non-renormalizability, the model predictions are usually compromising with the 
regularization method employed. The regularization cannot be removed from the results 
and consequently, the physical implications are crucially dependent on the adopted regular- 



35 



ization technique. It is usual practice consider the regularization as part of the model - see 
for example j|5, O, Is]. However, depending on the adopted regularization method, physical 
amplitudes may emerge from the calculations ambiguous and symmetry violating. Unde- 
fined quantities arising from divergences are fixed in the parametrization of the model by 
adjustments in regularization parameters. In veiw of this, one has the potential problem that 
the results of calculations are not real predictions, but particular choices for the involved 
arbitrariness (ambiguities) H] that lead to results that might change when different 
choices are made. 

Our investigation in the present paper focused on avoiding as much as possible explicit 
evaluations of divergent quantities. Our discussion consisted in basically two steps. In the 
first step we obtained physical amplitudes free from ambiguities and symmetry preserving. 
Integrands of divergent Feynman integrals were manipulated such that all the dependence on 
internal (arbitrary) momenta is left in terms that lead to finite integrals, which are then inte- 
grated free from the regularizations effects, while purely divergent objects (combinations of 
divergent integrals) can be clearly identified. Invoking symmetry constraints and demanding 
elimination of ambiguities lead to what we called consistency relations (CR), which require 
definite values for the divergent objects. All one- and two-point functions of the model, 
in which the highest degrees of divergence reside, emerge from the calculations free from 
ambiguities and respecting symmetry constraints. In a second step, a parametrization was 
introduced to eliminate divergences. We considered a parametrization of the model where 
the remaining divergent integrals I qua d (M 2 ) and I\ og (M 2 ), need to be eliminated by fixing 
phenomenological quantities. We used general scale properties of I qua d(M 2 ) and I\ og (M 2 ) 
in order to make explicit the freedom one always has when choosing a specific regularization 
through the constant parameter Ci- All this was done without explicit evaluation of the 
divergent integrals. 

Of course, we could also have chosen a specific regularization for evaluating explicitly 
Iquad (M 2 ) and /log (M 2 ) and put our calculation in close connection with the traditional 
regularization methods. As a result, both I qua d (M 2 ) and I\ og (M 2 ) become a function of a 
regularization parameter A and would also depend on the form chosen for the regularization 
function G(k 2 /A 2 ). In this way, for different regularizations one would have different values 
for C*2- We have shown that without making any choice one can fix this C<i by simply choosing 
its critical value. An important point to be noted, we reiterate, is that all manipulations 
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done prior to such a choice have been made guarantee amplitudes that are free from loop 
momenta ambiguities and that preserve the symmetries of the model. 

The assertion that the formulation of the regularization of the NJL model presented 
here is predictive is to be understood in the sense that no arbitrary choices were made in 
intermediate steps up the stage of calculating phenomenological quantities, like the pion 
and sigma masses and coupling constants. Even at this last stage, the existence of a critical 
value for C2 that leads to good values for phenomenological quantities is gratifying. 

There is a variety of applications that can be made using the general formulation presented 



here. Like the recent application to clustered hadronic matter 32J and color superconduc 



tivity in high density quark matter 33j, we envisage great potential for the study of heavy 



flavor in quark matter. The SU(3) version of NJL model constitutes a natural candidate to 
the application of the formulation presented here. Work along this line is presently under 
way [2^. 
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